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Abstract 

A theorem on the existence of exactly hmit cycles around a critical 
point for the Lienard system x + f (x) x + g (x) = is proved. An 
alogrithm on the determination of a desired number of limit cycles for 
this system has been considered which might become relevant for a 
Lienard system with incomplete data. 
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1 Introduction 



Limit cycles are isolated closed curves in an autonomous system in a phase 
plane. Determination of shape and number of limit cycles has been a challeng- 



ing problem in the theory of autonomous systems. Lienard system has been a 
field of active interest in recent past because of its relevance in various physical 
and mathematical problem [[5] — [S]]. Recently non-smooth Lienard systems 
even allowing discontinuities are also being studied. 
Here we consider the Lienard equation of the type 

X + f (x) i + g (x) = 0. (1) 

The Lienard equation ([T]) can be written as a non-standard autonomous system 

x = y- F{x), y = -g{x) (2) 

where F (x) = f (u) du. The phase plane defined by is known as Lien- 
ard plane. Lienard gave a uniqueness theorem [[3], [S]] for periodic cycles for a 
general class of equations when F (x) is an odd function and satisfies a mono- 
tonicity condition as x — oo. A challenging problem is the determination of 
the number of limit cycles for a given poljTiomial F (x) of degree (m) for the 
system ([2]) [[2], [5], [8]]. Recently we have presented a new method for proving 
the existence of exactly two limit cycles of a Lienard system Recall that 
the proof of Lienard theorem depends on the existence of an odd function F (x) 
with zeros at x = and x = ±a (a > 0) and that F (x) > for x > a and 
tends to oo as X — )■ oo. To weaker this assumption, we note at first that the 
existence of a limit cycle is still assured if there exists a value a > a (called an 
efficient upper estimate of the amplitude of the limit cycle) such that F (x) is 
increasing for a < x < a < Li, where Li is the first extremum of F (x) , x > a. 
Based on this observation we are then able to generalize the standard theorem 
for the existence of exactly two limit cycles. Our theorem not only extends the 
class of F (x) considered by Odani [[6], [7]], but also that of the more recent 
work of Chen et al [H] [See P] for more details]. 

In the present paper we prove the theorem for the existence of exactly 
limit cycles for the system ([1]). In the second part of the paper we present 
an algorithm to generate any desired number of limit cycles around the ori- 
gin, which is the only critical point for the system ([1]). Limit cycles represent 
an important class of nonlinear periodic oscillations. Existence of such non- 
linear periodic cycles have been established in various natural and biological 
systems [[3], [9], [1]]. It is well known that mammalian heartbeats may follow 
a non-linear oscillatory patterns under certain (physiological) constraints [1]. 
However, sometimes it becomes very difficult to obtain total information about 
a nonlinear system due to various natural constraints, as a result of which we 
obtain only a partial or incomplete data [4j. Our objective is to fill up those 
gaps and construct a Lienard system that may be considered to model the 
dynamics of the missing part of the phenomena in an efficient manner. 

To state this in another way, let us suppose that the Lienard system is 
defined only on a bounded region [— ai,ai], ai > having one (or at most a 



finite number of) limit cycles in that region. Our aim is to develop an algorithm 
to extend the Lienard system minimally throughout the plane accommodating 
a given number of limit cycles in the extended region. By minimal extension 
we mean that the graph [x, F (x)), of the function F which is initially defined 
only in < ai is extended beyond the line x = ai iteratively as an action 
induced by two suitably chosen functions {x) and H (x) so that (j) acts on 
the abscissa x and H acts on the ordinate F [x) respectively. Accordingly the 
desired extension F (x) of F {x), x > ai is realized as H o F (x) = F o (p[x). 
The choice of (p and H is motivated by theorem [2] so that the extension F 
satisfies the conditions of the said theorem. It turns out that can simply 
be a bijective function, while H may be any monotonic function admitting 
a < L (c.f. equation (Q), L being the unique extremum of F (x), x G [ai, 02], 



The paper is organized as follows. In section [2] we introduced our notations. 
In section |3] we have proved an extension of the theorem in [2J for existence 
of exactly limit cycles in the Lienard equation. In section H] we present the 
construction by which we can get a system of the form ([T]) having any desired 
number of limit cycles around a single critical point. Examples in support of 
this algorithm are studied in section O 



We recall that [3] by symmetry of paths, a typical phase path YQY' of the 
system becomes a limit cycle iff OY = OY'. 
We consider. 



F(ai) = F(a2) =0. 



2 Notations 





and 





YQY' 



It follows that 



dv = ydy + gdx = Fdy 



(5) 



so that 




(6) 



We define. 




g (u) du 



y+ (o) = oy, 



y. (0) = OY' 



Figure 1: Typical path for the Lienard theorem 



and let Q has coordinates (a, F (a)). Let a' and a" be respectively two positive 
roots of the equations 

G{a) = \yl{Q)-\F\a) (7) 
and G (a) = \yl (0) - (a) (8) 

Also let, 

a = max {a', a"} (9) 

We show that Vyqy' has a simple zero at an a < a [2] for the system (j2]) 
in Lienard theorem. It turns out that a provides an efficient estimate of the 
amplitude of the unique limit cycle of the Van der Pol equation [2J. This result 
has been extended in |2j for the existence of exactly two limit cycles as stated 
in the following theorem. 

Theorem 1 Let f and g be two functions satisfying the following 
properties. 

(i) f and g are continuous; 

(ii) F and g are odd functions and g (x) > for x > 0.; 
(Hi) F has +ve simple zeros only at x = ai, x = a2 for some ai > and 
some a2 > a, a being defined by (0) and a < L, where L is the first 
local maxima of F (x) in [ai , 02] ; 
(iv) F is monotonia increasing in ai < x < a and F {x) — ?■ —00 as x ^ 00 
monotonically for x > 02] 
Then the equation ([Ij) has exactly two limit cycles around the origin. 

It has been shown [2] that these two limit cycles are simple in the sense that 
neither can bifurcate under any small perturbation satisfying the conditions 



of theorem [H The existence of a satisfying an equation of the form ([9]) ensures 
the existence of two distinct hmit cycles. 



We generahze theorem [T] as follows. 

Theorem 2 Let f and g be two functions satisfying the following 
properties. 

(z) / and g are continuous; 

{ii) F and g are odd functions and g{x) > Q for x > 0.; 
{Hi) F has N number of +ve simple zeros only atx = ai,i = l,2,...,N 
where < ai < a2 < ■ ■ ■ < such that in each interval Ii = [oj, aj+i], 
i = 1,2, . . . , N — 1, there exists cxi, satisfying properties given by 
such that Ui < Lj where Li is the unique extremum in Ii, 
i = 1, . . . , N — 2 and L^v-i, the first local extremum in [a^^i, ajy]. 
(iv) F is monotonic in ai < x < (Xi^l i and |F (x)| — )■ oo as x — )• oo 
monotonically for x > Un- 
Then the equation ([I]) has exactly N limit cycles around the origin, all are sim- 
ple. 



Proof. We shall prove the theorem by showing the result that each limit 
cycle intersects the x— axis at a point lying in the open interval 



3 Existence of Exactly limit cycles for 
Lienard System 



y 




Figure 2: 



i = 0,1,2, . . . , N — 1, where ctQ = Lq is the local minima of F (x) in [0, ai]. 
By Lienard theorem and theorem [T] it follows that the result is true for = 1 
and N = 2. We shall now prove the theorem by induction. We assume that 
the theorem is true for N = n — 1 and we shall prove that it is true for N = n. 
We prove the theorem by taking n as an odd +ve integer so that {n — 1) is 
even. The case for which n is even can similarly be proved and so is omitted. 
It can be shown that |3], Vyqy' changes its sign from +ve to —ve as Q moves 
out of Ai (ai,0) along the curve y = F (x) and hence vanishes there due to 
its continuity and generates the first limit cycle around the origin. Next, in 
[2] we see Vyqy' again changes its sign from —ve to +ve and generates the 
second limit cycle around the first. Also, we see that for existence of second 
limit cycle we need the existence of the point a, which we denote here as cvi. 

Since by induction hypothesis the theorem is true for = n — 1, so it 
follows that in each and every interval {ak, cxk+i], k = 0,1,2, ... ,n — 2 the 
system ([2]) has a limit cycle and the outermost limit cycle cuts the x— axis 
somewhere in {dn~i, oo). Also Vyqy' changes its sign alternately as the point 
Q moves out of a^'s, i = 1,2, . . . ,n — 1. Since {n — 1) is even, it follows that 
Vyqy' changes its sign from +ve to —ve as Q moves out of a„_2 along the 
curve y = F (x). Since there is only one limit cycle in the region oo), so 

it is clear that Vyqy' must change its sign from —ve to +ve once and only once 
as Q moves out of An-i (a„_i,0) along the curve y = F{x). Also it follows 
that once Vyqy' becomes +ve so that it does not vanish further, otherwise 
we would get one more limit cycle, contradicting the hypothesis so that total 
number of limit cycle become n. We now try to find an estimate of a for which 
Vyqy' vanishes for the last time. 

We shall now prove that the result is true for N = n and so we assume that 
all the hypotheses or conditions of this theorem are true for N = n. So, we get 
one more point a„ and another root a„, ensuring the fact that Vyqy' vanishes 
as Q moves out of v4„_i through the curve y = F (x), thus accommodating a 
unique limit cycle in the interval (q;„_i,q;„]. 

By the result discussed so far it follows that Vyqy' > when a lies in 
certain suitable small right neighbourhood of cxn-i- We shall prove that Vyqy' 
ultimately becomes —ve and remains —ve as Q moves out of An {an, 0) along 
the curve y = F (x) generating the unique limit cycle and hence proving the 
required result for iV = n. 

We draw straight line segments XkXj^, k = 1,2,3, ... ,n, passing through 
Ak and parallel to ?/-axis as shown in figure [2l For convenience, we shall call 
the points X„, X'n, Y, Y' as B, B',Xq, X'q respectively. We write the curves 
^k = Xk_iXk, T'j, = X'i^X'j^_^, k = 1,2,3, ... ,n 

so that 

n n n 

YQY' = XoQX', = ETfc + XnQX'^ + = E (r^ + n.) + BQB' 

k=l k=l k=l 



and 



Vyqy' = E (^r, + ^r' ) + ^bqe' 



(10) 



k=l 



We shall prove the result through the following steps. 



Step (A) : As Q moves out of along A„C, Vr^. + Vr'^ is +ve and 
monotonia decreasing for odd k. 

We choose two points Q{a,F{a)) and Qi{ai,F{ai)) on the curve of 
F (x), where ai > a > a„. Let YQY' and YiQiY( be two phase paths through 
Q and Qi respectively. We have already taken Y = Xq, Y' = Xq, B = Xn 
and B' = X'^. We now take Fi = Zq, = Z^, B^ = Z„, B[ = Z'^ and Z^Z'^ 
as the extension of the line segment X^X'^ V k. Also we write = 
and = A^. If k is odd, then on the segments and A^ we have y > 0, 

F (x) < and y - F (x) > 0. Now, 

0<[?/-i^(^)]r, <b-^(^)k- 
Since (z) > for x > so we have 



So, by ([2]) we get 



-9 \x) 



_y-F[x)_ 
dy 



dx 



< 



< 



-9 (x) 



y — F [x 
dy 



< 0. 



dx 



< 0. 



Therefore by flTT]) we have 



(11) 



Tfc Tfc Afc Afe 

Since F (x) and rfi/ = ydt = —g (x) dt are both —ve along A^ for odd k, so we 
have 



F dy>0. 



(12) 



Next, on the segments and A'^ we have y < 0, F (x) < and y — F{x) < 
0. Now, 

0> [2/-F(x)]r, >[y-Fix)]^, 

k k 

So, by ([2]) we get 



dy 

dx 



> 



dy 
dx 



> 0. 



(13) 



Therefore by f[T3|) we have 

Vr'^ = lFdy= I i-F) ^dx > J (-F) ^dx = j F dy = V, 



-A. 



Since F (x) and dy = ydt = —g (x) dt are both —ve along for odd k, so we 
have 



^r; > Va' = Fdy>0. 



(14) 



From (I12p and (I14p we have 



Vr, + Vr'^ > Fa, + ^a', > 0. 

Therefore Vp,. + Vr^ is +ve and monotone decreasing as the point Q moves out 
of A„ along A„C. 



Step (B) : As Q moves out from along A„C, Vp^ + Vrj, is —ve and 
monotonia increasing for even k. 

On the segments Tk and Ajt we have y > 0, F (x) > and y — F (x) > 0. 
Now, 

0<[2/-F(x)]r^<[y-F(x)]^^. 
Since (x) > for x > so we have 



So, by ([2]) we get 



-9 [X) 
J/ - F(x)Jp 

dy 



dx 



< 



< 



-9 [X) 
_y- F{x)_ 

dy 



< 0. 



A. 



dx 



< 0. 



(15) 



Therefore by (fT5|) we have 



F rfy = / F^rfx < /" F^dx = I Fdy = Va^. 



Tfe Tfc Afc Afe 

Since F (x) > and dy = ydt = —g (x) dt < along for even k, so we have 



Vr, < Va, 



F dy <0. 



(16) 



Next, on the segments and A'^ we have y < 0, F (x) > and y~F{x) < 
0. Now, 

0> [?/-F(x)U > b-F(x)]^, 



so that by ([2]) we get 



'dy 

dx 



> 



'dy_' 

dx 



> 0. 



(17) 



Therefore by (fT7|) we have 

Since F (x) > and dy = ydt = —g (x) dt < along for even k, so we have 

Vn < = [ Fdy<0. (18) 



K 

From fflGjl and (fTSjl we have 

^r, + ^r^ < Va, + Vj,'^ < 0. 

Therefore V^^. + Vp'^. is —ve and monotone increasing as the point Q moves out 
of An along AnC. 

Step (C) : Vbqb' is — ^^e and monotone decreasing and tends to — oo 
as Q tends to infinity along AnC. 

On BQB' and BiQiB[ we have F (x) > 0. We draw BHi and 
parallel to x-axis. Since F {x) > and dy = ydt = —g (x) dt < along 
BiQiB[ so 

Vb^q,b[= J Fdy< j Fdy. 

BiQiB[ HiQiH[ 

Since [F {x)]jj^q_^jj,^ > [F {x)]^q^, and F (x) > 0, dy = ydt = —g{x)dt < 
along BQB' so we have 

Vb,q,b[ < J Fdy< J Fdy = Vbqb' < 0. (19) 

HiQiH[ BQB' 

Let 5" be a point on y = F (x), to the right of An and let BQB' be an 
arbitrary path, with Q to the right of S. The straight line PSNP' is parallel 
to the y-axis. Since F (x) > and dy = ydt = —g (x) dt < along BQB' and 
PQP' is a part of BQB' so we have 

Vbqb'= J Fdy = - J Fdy<- J Fdy. (20) 

BQB' B'QB P'QP 

By hypothesis F is monotone increasing on AnC and so F (x) > NS on PQP' 
and hence (!20l) gives 

Vbqb' < - y NS-dy = -NS j dy = -NS ■ PP' < NS ■ NP. 

P'QP P'QP 

But as Q goes to infinity towards the right, so NP oo and hence by the 
above relation it follows that Vbqb' —oo. 



step (D) 



From steps (A) and (B) it follows that ^ (Vrf. + Vr' ) in flTU]) is bounded. 

Therefore as Q moves to infinity from the right of An ultimately the quantity 
Vbqb' dominates and hence Vyqy' monotonically decreases to — cxd to the right 
of An- The monotone decreasing nature of Vbqb' inherits the same nature to 
Vyqy' as Q moves out of An along the curve y = F {x). 

By the construction of a„ it is clear that Vyqy' > at a point on the left 
of An and ultimately it becomes —ve when the point Q is at the right of An- 
So, by monotonic decreasing nature of Vyqy' it can vanish only once as the 
point Q moves out of An along the curve y = F{x). Thus, there is a unique 
path for which Vyqy' = 0. By (|6]) and symmetry of the path it follows that 
the path is closed and the proof is complete. 

4 Construction of a Lienard System with 
Desired Number of Limit Cycles 

We now present an algorithm by which we can form a Lienard system with 
as many limit cycles as required. We present the technique for two limit cycles 
around a single critical point. This technique can similarly be extended for 
n number of limit cycles. As stated in the introduction, this algorithm is 
expected to become relevant in a physical model with partial or incomplete 
data. 

Suppose in a given physical or dynamical problem, the function F of the 
Lienard equation ([1]) is well defined only with a finite interval [— ai, ai] denot- 
ing F (x) = /i (x) for X e [— ai, Oi] and satisfying the conditions: 

(i) /i is a continuous odd function having only one +ve zero ai 

{ii) xfi (x) < W X E [— fli, ai] 

{in) fi has a unique local minimum at the point Lq within (oq, ai), Uq = 0. 
Suppose it is also known that the system has a limit cycle just outside the in- 
terval [— ai,ai]. We have no information of F (x) beyond the interval. Our 
aim is to develop an algorithm to determine a function /2 as a restriction of 
F in an interval of the form [01,02] so that it satisfies the conditions of the 
theorem [2], ensuring the second limit cycle outside [01,02]. Now to determine 
/2 precisely from the information of /i in [oo,Oi] we need to define two func- 
tions 01 and Hi so that we get the abscissa and ordinates of /2 in the interval 
[oi, 02] respectively. The choice of 0i is motivated by Odani's Choice function 
[6] (c.f. remarks [1] for further details). The functions 0i and Hi are defined as 
follows: 



the functions (pn and are bijective such that 

4>iL ■ [oo, ^o] [ai,Li] , (piL (ao) = ai, (piL {Lq) = Li 
(j)iR : [Lo, ai] -> [Li, 02] , (t)iR (Lq) = Li, (f)iR (oi) = 02 



and Hi is monotone decreasing on [0, /i (i^o)] such that 

^lo/i :=/2o0i. (21) 



To make the definition fl2Tl) exphcit we define at first two monotone functions 
/22, and /2j:j and then introduce Hi parametrically by the help of two monotone 
decreasing functions Hi^ and Hir on [0, fi (Lq)] as 

HiL ■■ fiL (s) -> /2L (s)) , s e [ao, Lo] 
HiR : /lij (s) ^ f^R (s)) , s G [Lo, ai] 

U U) ^ I HiL (x) , if a; = fiL (s) , s G [oq, Lq] 
^ ^ ' \ /fiij (x) , if X = fiR (s) , s G [Lq, ai] 

The choice of is made on the basis of fi (x) defined on [— ai,ai] and the 
second zero 02 of F (x) that must he close to but nevertheless, less than the 
expected amplitude of the second limit cycle. We define the functions /2L and 

f2R as 

/2L : <PiL (s) HiL ifiL (s)) , s G [ao, Lq] 
f2R ■ <PlR (s) HiR {fiR (s)) , s G [Lq, ai] . 

We should note that in the definition of (pn and (piR we have used the condi- 
tions (oq) = ai, 4>iL (Lq) = Li and (piR (Lq) = Li, (oi) = a2. We could 
also have used the conditions 0il (ao) = Li, (piL (Lq) = ai and (Lo) = 
a-2,4>iR{ci-i) = Li instead, but in that case the function Hi and H2 must be 
monotone increasing. 

If X G [ai, Li], then x = ^i^, (s) for some s G [ao, Lq]. Therefore, 

/2L (X) = f2L (01L (S)) = //iL (/iL (s)) = /2L (s)) = fL (^) • 

So, 

/2L = /2*L- 

Next, if X G [Li, a2], then x = (piR (s) for some s G [Lo, ai]. Therefore, 

f2R (X) = f2R {<PlR (S)) = HiR ifiR (S)) = {<PiR (S)) = (x) . 

So, 

f2R = fm- 



Thus the unknown functions /2L and f2R can be expressed by known func- 
tions and so that we have 



Next we construct the restriction /a of the function F in [02,03] having 
unique local maximum at L2 (say) in (02,03). We assume two bijective func- 
tions 

02L : K, Li] [02, L2] , 02L (ai) = a2, 4>2L (^1) = L2 
and 4)2R : [Li, 02] ^- [L2, 03] , ^2R (-f/i) = ^2, 02i? (02) = 03 

and two more functions f^i^ and /|^. We define two monotone decreasing 
functions H2L and H2R on [0, /2 (-^1)] parametrically as 




/2L (a:) , X e [oi, Li] 

f2R{x), xe[Li,a2] 

fkix), xe[ai,Li] 

fmix): xe[Li,a2] 



H2L : /2L (S) ^ {^2L (S)) 
H2R: f2R{s) ^ f;R{<P2R{s)) 



s e [oi,Li] 
s e [^1,02] . 



We define 



fsL ■■ (f>2L (s) H2L (/2L (s)) , 

Ur ■■ hR (s) H2R {f2R (s)) . 



s e [ai,Li] 

s e [Li,02] 



so that as shown above we have 



fsL = /ai and f^R = J*^. 



Therefore, 




fsR ix) , 



/3i? (■«) , 



X e [02,1/2] 
X e [^2,03] 

x e [02, L2] 
X e [^2,03] 



We observe that 



fsL : 02L (01L (s)) -> il2L (/2L (01L (s))) , S € [oq, Lq] 

/si? : 02i? (01K (s)) ^- i/2R (/2R (^^'iR (s))) , s e [Lo, Oi] . 



We can similarly proceed and construct all the restrictions fk of the function 
F in [ofe_i, Ofe] for /c = 4, 5, 6, . . . , A*" so that the corresponding Lienard system 



have exactly limit cycles. Thus an incomplete Lienard system can be ex- 
tended iteratively over larger and larger intervals of x, having as many (simple) 
limit cycles as desired. We note, however, that the choice of iterated functions 
has as large arbitrariness except for the required minimal conditions of mono- 
tonicity satisfying theorem [2l The number of limit cycles for each such choices 
remain invariant. The problem of reconstructing data with a given number 
of limit cycles and having specified shapes is left for future study. We now 
illustrate the above construction by the following examples. 

5 Examples 

We now present some examples following the construction described in 
section HJ Here the figures has been drawn using Mathematica 5.1. 

Example 1 Let ai = 0.2, a2 = 0.5 and 



f\ {x) = 0.15 - 0.25^/ 1 - ^ ^ , -0.2 <x< 0.2. 

Here, Lq = 0.1. Let Li = 0.3. Let us choose 



2 



/2*^(x) = -0.15 + 0.25Wl 



/*^(a;) = -0.15 + 0.25Wl 



(x-0.3) 
(0.125)^ 

(a; -0.3)^ 
(0.25)^ 



Also, let 



\l {s) = ^As^ + B. 



To determine the unknown parameters A and B we assume that (pn (oq) = ai, 
(f)^j^ (Lo) = Li. Then A = 5 and B = 0.04. Next, let 



(f)iR (s) = VA's^ + B' 
with (j)i{i{Lo) = Li and 0ii?(ai) = 02- Then, A' = — and B' = . Then 



following the algorithm in section\^ we have 

f2L = fL [ai,Li] 
and f2R = f2R in [Li, 02] 

so that 

f2L (x) , X e [ai,Li] 



3 300 



f2{x) 



f2R (x) , X e [^1,02] 



We now define 



f2{x) 

4 



< X < ai 
ai < X < a2 



{x — 0.5) , X > a2 



Fix) 



to make -F+ continuously differentiable in [0, 00). The last part of the function 
F+ is taken to make monotone decreasing for x > 02 so that the function F 
defined below satisfy the condition that \F {x)\ — t- 00 as x 00 monotonically 
for X > 02- We take 

' F+{x), x>0 
F_{x) , X <Q 

We find two limit cycles which cross the +ve y-axis at the points (0, 0.26731065) 
and (0,0.5749823) respectively. So, y+{0) = y_ (0) = 0.26731065 and ai = 
0.254219124. So, the conditions ai < Li are satisfied in this example. Thus 
the existence of limit cycles are ensured by the theorem\E with g{x) = x estab- 
lishing the construction in section ^ The limit cycles alongwith the curve of 
F [x) has been shown in figure 




Figure 3: Limit cycles for the system in Example [T] with the curve of F {x) 



Remark 1 From condition (C2) in \B\ we see that 

g{(j)k (x))(/)^ (x) >g{x). 

If g (x) = X, then it gives 

(pk (x) (x) >x. 



Thus, in example U\ if we take (pk (s) = 0il (s) = a/As^ + B, then the above 
inequality gives 

nr^ 2As 

y/As^ + B- , > s 

2yjAs^ + B - 

i.e., As> s 

i.e., A>1. 

By the definition of f2i and f2 it follows that the remaining part of the condi- 
tion (C2) is satisfied if 

and in particular 

|F(0i^(O))| = |F(O)|. 
Since, (f)iL (s) G [oi, 02] and s G [ao, Lq] so it gives 

|/2(0il(s))| > 1/1 sG [ao,Lo] 

^■e., > se[a^,L^]. (22) 

Next, in particular the equality occurs at s = Oq = and so we have 

(0))| = |F(0)| 
=^ |//il(/il(0))|=0 
=^ |i/i^(0)| =0 

=^ H,L (0) = (23) 
since F (0) = and fiL (0) = 0. By our construction we also see 

s- Hi{.s) <0 Vs. 

Thus, (f)iL behaves like choice function described by Odani. Here, the condition 
{\22i) does not hold for the system discussed in example [H In fact, here 

\HiL{fiL{s))\ < \f^L{s)\ se[ao,Lo]. 

However, the conditions {viz. cti < Li, etc.) of theorem\^ are satisfied ensuring 
the existence of exactly two limit cycles. This shows that the theorem\^ and the 
construction presented above covers a larger class of functions F than those 
covered in [6\. The equality in (\2'3i occurs in example U\ only at the point 
s = ao = 0. However, the equality can occur at points where s ^ a^. We 
present example below to show this kind of behaviour. 



Remark 2 The function f2 in exampleUlis obtained from fi by reflection 
and translation along x-axis. However, it is clear from the construction of 
section\^ there is a plenty of freedom in the possible extensions of fi having 



a fixed number of limit cycles, as illustrated in examples IB md O In these 
examples we consider more general transformations so that the limit cycles 
are obtained having amplitudes close to those expected from the given physical 
{dynamical} problem. 

Example 2 In this example our target is to construct an example in which 



' I (x - 0.144)^ 

0.055518 - O.OSa/ 1 - ^ „ „ . , < x < 0.144 



0.04 



0.148506 - 0. 172988 Jl 



0.0910146 + 0.0209854 J 1 



(x -0.144)^ 



(0.206686)^ 



, 0.144 <x< 0.34 



(x - 0.407)' 

(0.06751554) 

(x - 0.407)^ 

(0.125376)^ 
^ -3.0000372 (x- 0.5), x > 0.5. 



2 ' 



0.34 < X < 0.407 



-0.2280727 + 0.340073 Wl 



, 0.407 < X < 0.5 



and 



F{x) 



F+ (x) , X > 
F_ (x) , X < 



Here, ai = 0.2, 02 = 0.5, Lq = 0.144 and Li = 0.407. It is easy to show that 
(piL (s) = V4.974392361- s2 + 0.0625 and (piR (s) = V2.019706 ■ s2 + 0.12376838. 

Here, 



fi {x) = < 



(x- 0.144)^ 

0.055518 - 0.08-1/1-^ „ „. , < x < 0.144 



0.04 



0.148506 - 0.172988. 1 



(x - 0.144)' 



(0.206686) 



2 ' 



0.148506 - 0.172988. 1 



(x - 0.144)' 



(0.206686) 



2 ' 



0.144 < X < 0.2 



0.2 < X < 0.34 



/2 (X) = < 



0.0910146 + 0.0209854 J 1 



(x - 0.407)^ 
(0.06751554)' 

I 

\2 ' 



r, 0.34 < X < 0.407 



/ (x — 407r 

-0.2280727 + 0.340073^ / 1 - ^ '■ V, 0.407 < x < 0.5 

(0.125376)^ 



The second part of the condition (C2) m [6\ i.e., the condition (f^ does not 
hold. In fact, 

|^^il(/il(s))| < |/il(s)| in (0,0.05290111) 

and \H^L{fiL{s))\ > |/ii(s)| m (0.05290111,0.144). 



The equality occurs at s = and s = 0.05290111. Here, we get two limit 
cycles crossing the +ve y-axis at the points (0,0.29039755) and (0,0.567249) 
respectively so that y+ (0) = (0) = 0.29039755 and = 0.2892792083. 
Consequently, ai < Li and the other conditions of theorem IE with g (x) = x 
are satisfied justifying the existence of exactly two limit cycles. These two limit 
cycles alongwith the curve of F (x) has been shown in figure^ 
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Figure 4: Limit cycles for the system in Example [2] with the curve of F (x) . 



Example 3 We now consider an example involving three limit cycles by 
taking Oi = 0.1, 02 = 0.2, = 0.4 and 



(x — 051^ 

/i (x) = 0.04422166 - O.O8WI - ^ -0.1 <x< 0.1. 

^ ^ ^' (0.06)' - - 

Here Lq = 0.05 and let Li = 0.15, L2 = 0.3. We take 



f*^ (x) = -0.04422166 + 0.08y 1 

f*^ (x) = -0.04422166 + 0.08^ 

It is easy to construct 



(x-0.15)' 



(0.06)' 

(s-0.15)' 
(0.06)' 



[s) = V5s2 + 0.01, se[ao,Lo] 



Next, we take 



3L 



0.0043819183 - 0.03W1 



{x-0.3y 



(x) = 0.0043819183 - 0.03W 1 



(0.101084111) 

{x-0.3f 
(0.101084111) 



We can similarly construct 



hhis) = V^ = 2s, se[ai,Li] 
(s) = 2s, s e [Li, 02] . 



so that (j)2L (ai) = a2, (j)2L (Li) = L2, (f)2R {Li] 
define 



L2 and (t)2R{a2) = 03. We 



0.04422166 - 0.08 Wl 



(x-0.05)' 



(0.06)^ 



-0.04422166 + 0.08W 1 - 



(x-O.lS)' 



(0.06)^ 



(x-O.S)^ 



0.0043819183 - 0.03W 1 - 

(0.101084111) 
^ 2.0100758 (x - 0.4) , 0.4 < x. 



< X < 0.1 
0.1 < X < 0.2 

0.2 < X < 0.4 



2 ' 



Fix) 



F+{x) 
F_ fx) 



X > 
X < 



to make F continuously differentiable. We can easily calculate that ai = 
0.12418214965 and 03 = 0.2354818163 and consequently ai < U for i = 1,2. 
All the other conditions of theorem\^ with g (x) = x are satisfied and hence we 
get three distinct limit cycles as shown in figure\^ alongwith the curve of F {x) 
defined above. 

Remark 3 Here the function F is defined in such a manner that \ F (Lq) \ > 
|F(L2)| implying that /32 mentioned in Theorem 3 of \8l or in theorem 7.12, 
chapter 4 of the book does not exist and hence these theorems are not ap- 
plicable for the corresponding Lienard system. 
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